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abstract 

fhis  paper  studies  the  behavior  ol  the  (PARCORf 
(pARtial  CORrelaiion)  coefficients  of  the  stochastic 
gradient  adaptive  lattice  filter  in  response  to  a  comples 
linear  chirp  FM  signal  in  white  Gaussian  noise  A 
single-stage  model  for  the  behavior  of  the  coefficients 
IS  developed,  and  the  recovery  error  is  derived  Also  an 
accurate  model  for  a  two-stage  filter  is  derived 


I  INTRODUCTION 

The  performance  of  adaptive  filters  for  sinusoids  m 
white  Gaussian  no:sc  has  been  studied  in  (lt.j2|  for  the 
LMS  FIR  filter  and  [31  for  the  Lattice  Filter  The  irack-mg 
characteristics  for  a  complex  chirp  signal  have  been 
investigated  in  H|.|5l.[61.{7)  for  the  1-MS  and  RLS  Filters 
In  this  paper,  we  will  study  the  behavior  of  the  P.ARCOR 
(PARtial  CORrelation)  coefficients  and  the  recovery  error 
of  the  SG  (Stochastic  Gradient)  Adaptive  Lattice  Filter  m 
response  to  a  complex  linear  chirp  signal  in  white  Gaussian 
x  noise.  The  expected  values  of  the  optimal  PaRCOR 
coefficients  and  a  first-order,  single-stage  analytical  model 
of  the  expected  values  of  the  PaRCOR  coefficients  based  on 
the  SG  update  algorithm  arc  derived  and  compared  with 
simulations.  It  is  noted  that  the  optimal  model  contains  only 
the  chirp  sinusoid  component  at  iteration  k  However 
simulation  shows  that  both  the  chirp  sinusoid  component  at 
Iteration  k  and  at  iteration  1  are  present  This  is  the 
"shadow"  effect  observed  in  (8|.  The  single-stage  analytical 
model  will  be  used  to  explain  how  the  update  algorithm  of 
the  SG  Lattice  retains  this  "shadow"  component  A  simple 
single-stage  model  of  the  recovery  error,  which  measures 
the  ability  of  the  filter  to  extract  the  signal  from  the  input 
will  also  be  presented. 

A  detailed  two-stage  model  is  derived  for  the  cases  of 
noiseless  and  noisy  input  signals.  It  was  shown  in  15|  that 
even  though  the  RLS  algorithm  e.xhibiis  superior 
convergence  characteristics  over  the  LMS  algoninm.  ns 
tracking  performance  is  in  fact  inferior  to  that  of  LMS  This 
IS  because  convergence  and  tracking  are  different 
phenomena.  The  results  derived  in  this  paper  show  that  this 
IS  also  the  case  for  the  SG  Adaptive  1-attice  Filter,  where  the 
convergence  and  tracking  rates  are  defined  by  independent 
parameters.  Furthermore,  the  analysis  shows  that  the 
PARCOR  coefficients  are  reduced  from  their  optimal  values 


when  either  the  chirp  laic  or  the  input  nOise  is  uKteaved 
The  second  stage  PAKCOR  coefficient  is  dirocd  assuming 
that  the  first  coefficient  has  reached  iis  steady  stale  !he 
accuracy  of  ilie  model  is  verified  bv  computer  simuiaiHvn 

2  niL  STOCHA.ST1C  GRADIENT  LATriCF  TH.TFR 

The  lattice  filter  structure  to  be  considered  is  shown  m 
Figure  I  The  laiiicc  order  recursion  equations  are  given  bv 


Ci(k|n)  =  e,(k!n  -  I)  KUhCk  -  l(n  U  (D 

Cpfktn)  V  c^(k  -  tin  1)  -  Kj;  e,(kin  -  1 ,)  C-i 


Figure  1  l-atticc  Filter  Structure 
The  of  K mmirni^es  'be  mean  squared 
forward  prediction  error  E[cf{k|n)|  is  gi'en  by 

F:|e,(kln-  DcUk  -  l|n  -  Mi 

E(c;(k  -  l|n  -  1)1  (  t) 

where  we  have  assumed  that  E(|ci(k|n)i'|  -  E|lCs(k!n)|- j 
and  that  K'„(k)  ^  K!;  (k) 

To  evaluate  (it)  we  determine  the  forward  and 
backward  pred.cnon  errors  of  an  n-iap  transversal  filter 
This  is  based  on  the  assumption  that  the  expected  forward 
and  backward  prediction  errors  of  the  lattice  filter  at  stage  k 
are  statistically  equal  to  those  of  a  k -step  irans' ersal  filter 

3  OPTIMAL  PARCOR  COEFFICIENT 


The  sicnal  model  used  here  iv  that  of  a  complex  linear 
chirp  f  M  m  white  Gair  sian  noixe  The  input  signal  xfkl  ^ 

given  bv  /s  ^ 

x(k)  .  s(k)  .  n(k) 

where  n(k)  is  a  rero  mean  white  (laussian  noise  signal  and/^^^/y^ 
s(k)  IS  a  sinusoidal  linear  FM  signal  defined  as 
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'.\here  ui  IS  Ihe  base  Ife()uc(HH  'P  ts  Ihe  chirp  rate  <t>  ii  a 
random  phase  ihat  is  assumccj  lo  be  uniformly  disiribuied 
beiueen  -n  and  n.  P,  is  ihe  si^nai  pow,ei,  and  P„  is  the  noise 


power 

The  characiensncs  of  the  SG  Ijiiice  filtei  can  f>e 
deroed  bs  assummc  that  the  backward  and  forward 
prediction  errors  at  stage  n  are  statistically  equal  to  the 
respecioe  backwanl  and  forward  prediction  errors  of  an 
n-siaee  linear  predictor  The  optimal  PARCOK  coell<c>eni 
can  be  obtained  by  assuming  that  the  backward  and  forward 
errors  have  anained  their  optimal  values  This  can  be  shown 
(Figure  2)  to  be 


K;;‘’‘(ki  =  a*b 


a 


I  •  no 


i  fo 


where  a  v  e'*"  b  --  e'’‘ 


and  p 


P 

P„ 


The  update  algorithm  of  the  SG  lattice  algorithm  does 
not  permit  the  reflection  coefficients  to  anain  their  optimal 
values  for  the  non-siaiionary  chirp  FM  signal  Instead,  the 
shadowing"  effect  is  observed  Simulations  with  complex 
FM  signals  reveal  that  the  actual  reflection  coefficients  at 
each  Iteration  contain  two  frequencies  -  the  fundamental 
frequency  at  the  start  of  the  S'veep  and  the  swept  frequency 
at  that  iteration  The  next  section  will  derive  a  first-order 
approximauon  of  the  reflection  coefficient  characteristics 


MTHiUAi  i/r-n  r»fry(>**K'r 

Figure  2  Transfer  Function  of  the  Expected  Value  of  the 
Optimal  PaRCOR  Coefficients 

4.  SINGLE  STAGE  PARCOR  COEFFICIENT 


The  reflection  coefficient  update  equation  is  given  hy 
K„(k*  1)=  |1  -/5l)er,(k-  l|n-  l)|'|K„(k) 

♦  /Se,(kln  -  UCbfk  -  Ijn  -  I)  (ty 


where  (3  is  the  forgetting  factor  (equivalent  to  p  m  the  L.MS 
algorithm)  For  convenience  we  shall  drop  the  subscripts  m 
n  and  n- 1  Taking  expectation  on  both  sides,  the  neraied 
solution  of  (7)  IS  given  by 

t 

F|K(kv  I)|  =  K(0)Elf](l  -  l)|‘! 

t-  \ 


^  e,(k  -  r)e^(k  -  r  - 

f 


n[]ll  -/3|e^(k  -  s)l- 


iM 


Assuming  itiji  KiUj  -  0.  and  ncgtcviing  iffisw, 
highci  (8)  vail  t,K  evaluated  to  tve 


FiKtk 


fibP,  a  a'*' 
I  .  I'd  i  a 


Due  to  the  vummation  m  (8)  wliich  luio',  out  i<i  tK  a 
eromeuK.  venes  unlv  ihe  lundamenial  lenn  a  and  it,-- 
tetm  at  the  krlsi  iieiaiwin  a"'  irmain  Sins  vrivrv  ; 
vonbim  the  observed  Siinulation  lesuits  jv  n'. 

T  leoicv  ‘  ariil  J 


Figure  3  Transfer  Function  of  the  Expected  PAR(  OR 
Coefficients  from  Simulation 


Figure  4  Transfer  Function  of  the  Expected  P-XRf  OR 
Coefficients  from  Single-Stage  Model 

5.  SINGLE  STAGE  RECOVERY  ERROR 

The  recovery  error  m  this  Situation  mcj'urev  ihc 
amount  of  residual  signal  left  m  the  forward  error  XSe  shall 
derive  a  simplistic  single  stage  recovery  error  expression 
here  The  recovery  error  is  defined  as 

7(k|n)  =  edkln)  -  n(kj  i  tOi 

Defining 

r7a(k|r.)  =  e,(k|n  -  1 1  -  K;!(k)e^(k  -  l|n  -  1 1  -  niki  i  n  > 
we  can  rearrange  OO)  as 

Vf'vln)  -  >/o(F|n|  =  -  I  ,>(k)e^(k  -  !jn  -  1  !  O'* 

where 

I  „(k|  -  K„tk)  -  K2(k| 

We  wish  to  determine  the  quantuy  f.||i;(kin|  [ 

following  the  methodology  and  assumpnonv  of  ioj  ibt’ 
above  quaniitv  can  he  found  to  be 


,  1870 


E|lr?(kln>  -  'Mtvin)!*'! 


where 


M  = 


^  -  I  ■  _£'P„ _ 

M  -  a  i !  «•  pn|( !  ♦  ()(n  -1)1 

0Pi  I  »  p{n  »  1 ) 

2  1  ♦  p(n  -  f) 

l>  »pn|  -0 


[IJ) 


I  ♦  {){n  -  1) 

6  T^VO-STACE  ANALVSiS  (NOtSELESS) 


This  section  presents  a  two-stage  analysis  of  the  SG 
lattice  filter  assuming  that  the  input  contains  no  noise  The 
next  section  wii)  analyze  the  case  with  input  noise  For  a 
iHO-staee  model,  the  input  signal  is  given  by 

e,{k|0)  =  e^klO)  =  (14) 

Letting  a  =  e'***  and  b  =  e*".  we  get 

e,(k-riO)c;(k-r- 1|0)  =  P,ba‘-'/-a'<‘-'’  (IS) 
The  reflection  coefficient  update  formula  is  given  by 
K,{k*l)  =  (1  - /3T»,lK.,('f)  »  ^e,(k|0)e|.(k  -  t|0)(J6) 

Let  q  =  I  -  /fP <  Then  from  (15).  K,(k  *■  I)  can  be  evaluated 
to  be 

,  (11*11 

K,(k  V  1)  =/3P,b’a'‘a'-'^-^ - : - 9-^ 

a  -  q  (17) 

The  plot  of  K,  verses  iteration  number  is  shown  in  Figure  7. 
The  simulation  and  analytical  model  results  are  identical  in 
this  case 

From  (17)  we  see  (hat  K,(k*  I)  has  a  transient  and 
steady  state  component.  The  transient  component  is 

'  dependent  on  l  -  ^P, .  This  gives  a  convergence  time 
constant  of 

I 

°  /?P,  (18) 

^P» 

When  the  chirp  rate  is  zero,  the  term  a'  -  q  approaches 
unity  This  gives  Ki(k+  !)  its  maximum  magnitude.  With 
increasing  chirp  rate  this  term  reduces  K,(k*I)  in  the 

manner  shown  in  Figure  5,  where  a'  -  q  is  graphed  as  a 
function  of  chirp  rate  for  different  P's. 

To  derive  the  second  reflection  coefficient  Kj(k  *  1). 
we  shall  approximate  K|(k  +  1)  by  its  steady-state  value 

K,(kv  1)  H  )SP,b'a  ‘a-’«A^  (19) 

a  -  q 


Figure  5  f'loi  of  a  q  foi  a«!!c(fin  (j  and  chirp  lairc 

The  forward  and  the  backward  predicnon  cnon  of  ibr  fuM 
Stage  are  derived  from  the  /eroih  stage  by  ihe  follow  me 
Iterative  equations 

ei(ktl)  -  e,{klO)  K,(k>e^(k  -  HO) 

Cb(k|I)  ■*  e,(k  IjOl  K;  c,(k|0)  i:ui 

Defining  l  -  ^|c,(k|l)j’  .  q,  the  second  rellecnon 
coefficient  can  be  wnnen  as 

t  2 

Kjfk  ♦  1)  *  )9^c,(k  ■  rj!)eh{k  -  r  ~  l|I)q, 

..0  (21) 
Using  the  iteration  of  (20).  (2t)  can  be  evaluated  to  be 
K,(k  V  1)  =  I, 

where 

I,  =  /3P,b^9- 


Ij  =  I, 


_  ^  a'“ 


a-q 


q‘  'a 


a'^-  q, 


u 


a-q 


-  a'^  -  qV'a 


1,  i,  ,  1. 

(22) 

-  ' 

-  qi 

(-’.■'I 

-  ^ 

qa- 

qi 

(2*1) 

q‘-'3'‘*'  - 

q^-'a^ 

q‘a  ‘ 


qi 

q‘qa 


qa-q, 

q^^-'a'  -  qV'qa 


qa-q, 


q  -qi 

Examining  1,-1,  we  sec  that  the  convergence  of 
Kj(k+1)  is  dependent  both  on  q  and  q,  The  time  constant 
related  to  q  is  given  by  (18)  The  time  constant  related  to  q, 
IS  given  by 


/?P. 


1  - 


Ico'iy-  Zq  -  gp,  1 
1,^  -  l^ coiW  *  1 


(26) 


Figure  6  shows  a  plot  of  Kj  versus  iteration  number  Noie 
that  the  convergence  rate  of  the  second  stage  is  slower  than 
th'^i  of  the  first  The  steady-state  '■alue  of  Ki(k*l)  can  be 
derived  to  be 
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hm  {vCSt  +  1) 


a'-’-q 


a  -q  {a  -q)^ 


(27) 


(  .».r  |,..v|,l  r,,i  -  .fll.l-  .  «  0^  {M-«*  •  «  «#<“*• 


Figure  6  Real  Componeni  o/  the  Second  PARCOR 
Coefficient  (K2) 

7.  FIRST  STAGE  ANALYSIS  WITH  NOISE 


The  general  iterated  solution  for  the  first  reflection 
coefficient  is  given  by 
k-l 

K,(k+|)  =  eKk-r|0)eb(k-r- 1|0) 

■  nl  l-^e^{k-s|0)| 
(28) 


The  noise  terms  in  i  -  /3e^(k  -  s|0)  for  s=l  to  s=r-l  arc 
independent  of  each  other  and  of 

e,fk  -  r|0)e^  (k  -  r- l)|  1  - /3e|(k  - siO)|.  Thus 
E(Ki(k  -f  l)|can  be  wriacn  as 


k-l 


EiK,(k  +1)1  =  ^  X  E[eKk-r(0)eb(k-r  -  1|0) 

r-S 

j  l-^eb^k-r|0))]  n|l-^E|!cb(k-s|0)i')| 

(29) 


s=  1 


Letting  qn  =  (  1  -^Ele^(k-slO)|  =  I-^(P,  +  P„).thc 
summation  can  be  evaluated  to  be 


ElKifk  +  1)1  =  ^Pjb'a-’'- 


-'-gna' 

a’  -  qn 


l-m  -*•  2Pn) 
+  Pn) 


(30) 

We  see  that  the  reflection  coefficient  depends  upon  the 

signal  and  noise  by  the  factor  .i — and 

a  -  qn  1  -  ^(P,  +  Pn  ) 
the  convergence  constant  is  now  given  by 

1 

°  ^(P,  *P„)  (31) 

Figure  7  shows  the  plot  of  (30)  superimposed  on  ihe 
actual  simulation  -The  second  reflection  coefficient  iteration 


formula  has  also  been  denied  bui  shali  not  tie  ptrscr.ird 
due  to  space  considciations 
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Figurr  7  Real  Component  of  the  F.rsi  PaRCOH 
Coefficient  (Kl)  wmh  noise 


8  CONCLUSION 


In  this  paper,  we  have  studied  the  behavior  ol  the 
PARCOR  coefficients  of  the  SG  Lanice  Filter  for  a  comples 
linear  chirp  FM  in  white  noise  Independent  lerms  which 
define  the  tracking  and  convergence  rates  are  derived 
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